Interaction and heat exchange in 
two-component relativistic fluid composite 



Ernst Trojan and George V. Vlasov 
Moscow Institute of Physics and Technology 
PO Box 3, Moscow, 125080, Russia 

August 16, 2011 

Abstract 

A model of two-component relativistic fluid is considered, and the 
thermal nature of coupling between the fluid constituents is outlined. 
This thermal coupling is responsible for non-ideality of the fluid com- 
posite where the components are not fully independent. The inter- 
action between particles is reflected only in the equation of state of 
each component, but it deals nothing with the coupling between the 
fluid components and does not influence the hydrodynamic motion. 
A general form of two-fluid decomposition is formulated for arbitrary 
interacting system. 

1 Introduction 

Two-component and mult i- component relativistic fluid systems are often con- 
sidered in various problems of nuclear physics and astrophysics. As a rule, 
the matter is described in the frames of ideal fluid mechanics. However, the 
coupling between fluid components may give rise to non-ideality of the fluid 
composite. Relativistic hydrodynamics of non-ideal multi-component fluid 
p] involves complicated equations but it allows to operate with more precise 
theory. Particularly, the theory of two-component relativistic fluid has found 
important applications to hydrodynamics of relativistic superfluidity [2j [3] 
and nuclear hydrodynamics of proton- neutron systems [I]. The stress-energy 



tensor of a two-component fluid [5] 

T vu = n u fi a + s u <d a + P g ua (1) 

includes the particle number current n u and the flux of chemical potential 
/i CT corresponding to the first component, and the particle number current 
s v and the flux of chemical potentia Q a corresponding to the second compo- 
nent. In a relativistic superfluid the variables n u and [i a belong to the "cold" 
component, while the entropy density current s u and the temperature flux 
Go- belong to the "warm" component. 

The vectors in ([[]) are not independent but obey relationship 



V 



Q G J I 6" 



(2) 



with the coefficients 



1 dP 1 dP ldP 

F =-7T G = 7^ Q = ~lf ( 3 ) 

because the pressure P, in general, depends on three variables 

^ = -fffr e 2 = -e^ z 2 = -/j, u g u (4) 

Substituting (jBJ) in ([!]) we write the stress-energy tensor in the form 

T va = Fn„n a + GQ u Q a + 2Q^ u Q a + P g ua (5) 

The vanishing coefficient Q = will imply ideality of the two-fluid system 
whose stress-energy tensor (j5]) will be written in the form 

T va = Ffi u fi a + GQ u Q a + P g va (6) 

or 

T va = iinu u u a + Qsv u v a + P g va (7) 

It is a model of two ideal fluids without coupling. Velocities u v and v u can 
be different but there is no interference between them. 

So, non-zero coefficient Q implies that the motion of fluid components is 
not fully independent. The pressure P depends on the cross term z. This fact 
reflects dependence on the relative velocity between the fluid components 



u, = A /l-L- (8) 
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In the light of ([3]), it results in non-zero coefficient Q that implies coupling 
between the components and non-ideality of the two-fluid system: vectors n u 
and \i v are not collinear (as well as s v and Q u ). 

The coupling between the fluid components should not be mixed with 
the inter-particle interaction. The first component can be an ideal gas of 
free particles, the second component can be an ideal gas of free particles or 
thermal excitations (quasi-particles), but the fluid components are coupled as 
soon as the pressure P depends on the relative velocity w ([8]). It may look as 
an artifact of a plain mathematical trick without solid physical background: 
as soon as the relative flow w is introduced, the thermodynamical parameters 
will contain functional dependence on w. However, this idea was successfully 
developed in the Landau two-fluid model of superfluid helium and it was 
confirmed in experiments [6]. Although the source of coupling between the 
" cold" and " warm" components of relativistic superfluid is a consequence of 
interaction between the particles [3J, dependence on the relative velocity w 
is not evident in the field-theoretical approach to interacting many-particle 
system [21 [7j, and it had been a subject of dispute for a while [31 [8]. 

Researchers consider both variants - either a mixture of two ideal fluids, 
or a two-fluid composite where the coupling between components implies 
non-ideality of the whole fluid system. The first model (without coupling) 
is simple and it is often applied in astrophysics and nuclear physics. The 
latter model (with coupling) is much more complicated and it promises to 
give comprehensive description of hydrodynamic motion, although it is rarely 
used on account of its complexity. Of course, it is very desirable to take into 
account the role of coupling when we consider a strongly-interacting medium, 
like superfluid or nuclear matter. 

In the present paper we discuss the nature of coupling between the com- 
ponents of a two-fluid composite and emphasize that it arises from the heat 
exchange between them. The interaction between particles in a quantum 
many-body system will also result in non-ideality of the continuous medium. 
It is necessary to clarify what factor plays dominant role in the hydrody- 
namics of a two-component fluid: the interaction between particles or the 
thermal contact between fluid components? 

Standard relativistic units cu g ht = h = ks = I are used in the paper. 
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2 Thermal coupling between components 



Taking Equation (jSj), we immediately determine the energy density of a two- 
component fluid 

e = t ° = iWo + G0 o e o + 2g^ e + P g° (9) 

Non-ideality of the two-fluid composite, is due to non-zero term Q. In a 
relativistic superfluid this is introduced at the macroscopic level of hydrody- 
namics [5j. However, its appearance is a consequence of internal processes 
in superfluid system [9] and it implies dependence on the relative velocity 
between the "cold" and "warm" component w, that is confirmed in exper- 
imental research [6]. The similar term Q will appear in the stress-energy 
tensor of an arbitrary two-component fluid if its components are not fully in- 
dependent PQ. What process is responsible for non-ideality of this two-fluid 
composite? In order to understand how this non-ideality arises, let us turn 
from classical description to the quantum level. 

For an arbitrary classical fluid system with the stress-energy tensor (jSJ) 
and the energy density E (Q we can recognize the Hamiltonian H that 
determines the ground state |$) according to the stationary Schrodinger 
equation 

H\$) = E g m (10) 

with the ground-state energy 

E g = ($|#|$) = / T °d 3 r (11) 



Let us present the field variables \x v and Q u in the form [TD] 

fiu (x) = ^ Yl { e - ( fc ) Sfc e ~ ikx + < ( fc ) a I elkx ) ( 12 ) 

k 

(*) = 7 E (*) ik e ~ lkx + < (*) ^ elkx } ( 13 ) 

k 

with a unit polarization vector e u and form-factors X and Y, while the 
annihilation operators a k and b k satisfy commutation relations 

[a k , aj] = 5 kp b k ,bl = S kp (14) 
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and 



bk , bp 



a\ , dp 







(15) 



Substituting ffT21) . f[TH]) in ([2]), we find the Hamiltonian 



if = » <e = v £ (4a* + i) + ^ (i& + i) + ^ (aft + it, 

(16) 

where V is the volume of the system, F^, Gk and Qk are Fourier images of 
F (x), G (x) and Q (x), while O includes the terms with zero average in the 
light of (|T4|) - ffT5|) . Expression 



is no more than the operator of number of particles with energy 



and momentum k. Expression 



X 2 



N\ = b{b k 



is the operator of number of particles with energy 



e 



and momentum k. Quantity 



A 



k y 2 

2Qk 

XY 



(17) 



(18) 



(19) 



(20) 



(21) 



also plays the role of energy. It reveals the fact of coupling between the 
components because = is equivalent to Q = 0, and the Hamiltonian of 
this ideal system 

H° = E { £ * ( fi t fi * + + e * ( S & + I) } ( 22 ) 

corresponds to a two-component ideal fluid with the stress-energy tensor ([6]) 
or flU). 
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The Hamiltonian (j22p is already written in a diagonalized form. In or- 
der to diagonalize the Hamiltonian of non-ideal fluid (I16p let us apply the 
following transformation 



a k = cos 7] k a k + sin r] k f3 k 
b k = cos 9 k h + sin 6 k a k 



(23) 
(24) 



According to (JHJ) - (fT5]) . the operators a k and (3 k also satisfy commutation 
relations 



[a* , «J] 



5, 



kp 



[Oik , a p J 







(25) 
(26) 



Then, the Hamiltonian (1161) will be presented in a diagonalized form 
H = W 



^2 { ( £ t c °s 2 Vk + £ a k sin 2 9 k + 2£ k cos r] k sin 

k 

+ (e\ sin 2 r/ fc + e\ cos 2 + 2£ fc sin % cos fe ) ftj k } (27) 



where 



k 



(28) 



is the zero-point vibration energy, and parameters of transformation 
must satisfy condition 

2X k (cos rj k cos 9 k + sin r) k sin 6> fc ) + s s k sin (27/ fc ) + e£ sin (29 k ) = 

Particularly, at X k = condition ( 129]) is simplified 

e a k sin(2%) + e£ sin(20 fc ) = 

and its solution is sin T] k = sin 9 k = that implies no more than 

Sfc = otk h = (3 k 



(29) 
(30) 
(31) 



because the Hamiltonian ( 122]) is already given in a diagonalized form. 

Thus, condition (129]) imposed on transformation (1231)- (124]) determines a 
link between parameters t] k and However, it is not enough for their 
ultimate definition. An additional condition 



el = el (cos 2 i] k + sin 2 9 k ) = e k (sin 2 r] k + cos 2 9 k ) 
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(32) 



determines exact i] k and 9 k that allow to present the Hamiltonian (1271) in the 
following universal form 

H = Hideal + ^2 Uk _ PkPk) ( 33 ) 

k 

where 

H i( ieai = W + J2 (4 + A fe ) [a\a k + (34) 

k 

is the Hamiltonian of ideal system, while 

u k = A fc sin (7} k - 9 k ) (35) 

and 

A k = \ k sm(r] k + 9 k ) (36) 

The energy shift A k reflect no qualitative difference from the ideal system 
with Q = and X k = whose Hamiltonian is 

H = Wo + J2 4 (4 & k + PA) (37) 

k 

So, there is no qualitative difference between two systems with Hamiltonian 
f[3"4"|) and f[3"T|) . They both describe an ideal fluid. As soon as nonzero Q ^ 
appears in the stress-energy tensor (jSj), it implies non-ideality, and the 
relevant Hamiltonian (133]) includes additional nontrivial terms with a\a k — 
(3\f3 k . To understand their role let us appeal to the formalism of thermo field 
dynamics [H1II2]. Let us introduce the operation of tilde-conjugation 

a k = (3k (38) 

For operator 

h = ^ u k a\a k (39) 

k 



h = ^ w fe a[« fc = Wfe/3^/3 fc (40) 



its tilde-conjugated counterpart will be 
and we can rewrite the Hamiltonian (133]) in the following form 

H = H ld eal + h-h (41) 
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while the Hamiltonian H^eai is invariant under tilde-conjugation 

Hideal = W + ^2 ( £ k + (&\&k + ttl<5fc) = + ( £ k + ( 5 fc fi * + °i^k 



If we introduce thermal doublets 



(42) 



Ai 



OCk 



a. 



-at. 



(43) 



which satisfy commutation relations 



^4fc) A p 



'kp 



then, the Hamiltonian (j4T[) will be presented so 



or 



where 



and 



ideal 



+ 5^w fc ( .-I/, .4, + I 



ideal 



it 



1 



£*fa«j = + E K + A *) («*«* + ftp* 



(44) 



(45) 



(46) 



(47) 



(48) 



Again, the shift of the zero-point vibration energy (1481) reflects no qualitative 
difference between the Hamiltonians (l4"Tj) and (|37|) . However, the last term 
in the Hamiltonian ( 146]) makes essential difference from the ideal system with 
the Hamiltonian (1371) . 

We can take the Hamiltonian of non-ideal fluid in various forms, applying 
the tilde-conjugation (138|) and thermal doublets (|43l) . However, no transfor- 
mation can reduce this Hamiltonian to the form which corresponds to the 
Hamiltonian of an ideal fluid. The last term in ( |4"6]) is responsible for the 
coupling between the components. We can consider the Hamiltonian of non- 
ideal fluid in the form (|41|) . where the last two terms are responsible for the 
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coupling. In the Hamiltonian taken in the form (1331) the last terms, again, 
include the coupling. Now we can explain the very nature of the coupling 
between the components of non-ideal fluid. It is resulted from exchange of 
thermal excitations in the process of absorption of a quantum with the energy 
Uk and momentum k (action of operator or in the process of annihila- 
tion of a hole with the energy — and the same momentum k (action of 
operator ftk). The excitation and the hole with respect to the zero energy 
level of the ideal fluid E' Q ( 148|) . If the components are isolated and there 
is no heat exchange between them, then, Uk = and it is necessary to be 

h = h = that corresponds to an ideal fluid. Thus, the heat exchange is the 
source of coupling between the fluid components. The similar interpretation 
is known in hydrodynamics of superfluid helium where collective motion of 
thermal excitations is responsible for hydrodynamic behavior of the whole 
system [6]. Although this link between thermal excitations and macroscopic 
properties of the continuous medium is intuitively evident, we have proved 
it in the explicit analysis with formulas ( 1331) . (|41|) and (|46[) for the first time. 
However, it is still unclear what is the role of interaction between particles. 



3 Interaction between particles 



The pressure P and energy density E of a many-particle system are defined 
by standard formulas [T3] 

PV = -QlnZ 



and 



EV = & 



dlnZ 
dQ 



where the statistical sum Z is determined by formula 



InZ (e p ,fx, 6) = =F 



7 y 



(27T) 



d 3 p In 



1 ± exp 







(49) 
(50) 

(51) 



with the upper and low sign corresponding to fermions and bosons, and e p 
is the single-particle energy. For a system in motion at a velocity w the 
statistical sum is determined by formula [5] 



jV 



(2k) 



d 3 p In 



1 ± exp 



e 



(52) 
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and the distribution function is [2] 



f P (©) = — n (53) 

exp \e' p - p) /B ± 1 

where 

£p = £p-p-w (54) 

is the single-particle energy in the moving reference frame. 

Now consider a two-component fluid composite. When there is no relative 
flow between the components, the statistical sums of each component will be 
determined by the same formula (T5TT) in the co- moving reference frame, or 
by the same formula (1521) in the laboratory reference frame. What happens 
if relative motion between the fluid components occurs? Let us consider this 
two-fluid system in the reference frame co-moving the first component. Then, 
the statistical sum of the first component is defined by formula ( 1ST]) , that can 
be written briefly so 

Z a = Z(e;,pi a ,e) (55) 

The statistical sum of the second component is calculated, according to (|52|) . 
in the reference frame co-moving the first component: 

Z b = Z (e b p - p ■ w, fi b , 0) (56) 

where w implies the relative velocity between the first and second compo- 
nents. 

The total statistical sum of two-component system is 

Z = Z a Z b (57) 

Hence 

In Z = In (Z a Z b ) = \nZ a + In Z b (58) 

and, according to ( )49l) and ( |50|) . the pressure and energy density are additive 
quantities 

P = P a + Pb (59) 

E = E a + E b (60) 

How to complete a two-fluid decomposition for a system of interacting 
particles? It is convenient to apply the ideology of density functional the- 
ory that considers a system of free particles under the action of external 
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self-consistent field as a model of real interacting system. This ideology is 
developed in the nuclear mean-field approximation [TJ] where fluctuations 
of meson fields are neglected, and particles move independent in the mean 
fields, which themselves are generated self-consistently by the particles. This 
mean-field approximation allows to split the statistical sum 

Z = ZZ^ (61) 

into a product of the statistical sum of mean fields Z and statistical sum 
of particles Z^ where the latter is defined according to (15T]) . A set of self- 
consistent mean-field equations determine the effective mass M* and effective 
chemical potential //* which depend on the constant of interaction. As a 
result of P9|) . (150)) and (EE}, the pressure and energy density of the whole 
interacting system is 

P (e* p , //*, 0) = (e* p , //„ 6) + P (e* p , //„ 6) (62) 

E (e* p , /i„ 6) = (e* p , fi„ 6) + E (e* p , 6) (63) 

where P^ and E^ imply contributions of the particles, and P and E are 
contributions of the mean fields, while 

e* p = VP 2 + Ml (64) 

Again, following (14"§|) . (ISTJj) and ( 1o"5"j) - (l5"7]) . we determine the pressure and 
energy of the first component as 

P a = P (e* p a , fit, 6) = P^ (e* p , /i„ 6) + P (e* p , ^, 6) (65) 

E a = E (e; a , ft:, e) = ( e ;, ^, e) + e ( e ;, ^, e) (66) 

while the pressure and energy of the second component will be 

p b = p (e; b -p-w, & e) = p^ ( £ ; b - p- w, e)+P (e* p b - p- w, fil e) 

(67) 

E b = E (e* p b -p-w, fit, 6) = E* (e* p b -p-w, /it, Q)+E (e* p b -p-w, 6) 

(68) 

As soon as we know the thermodynamical functions of interacting system 
(l62|) - fl63|) . we can define the thermodynamical functions of the two compo- 
nents fl65l) - fl66l) and (J6"7]) - (J68"]) . thus making a two-fluid decomposition (J59|) - 

dnoD. 
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4 Conclusion 



Non-ideality of the two-component superfluid at the macroscopic level is ex- 
pressed in additional coefficient Q that appears in the stress-energy tensor 
(JSJ) and reflects its dependence on the relative flow w, resulting to coupling 
between the fluid components. The analysis in the frames of thermo field dy- 
namics has revealed that the Hamiltonian contains irreducible terms respon- 
sible for coupling and implying the heat exchange between the components. 
This effect of thermal contact is taken into account when the relative velocity 
between the components is introduced in the thermodynamical functions of 
a two-fluid system. 

One can consider the thermodynamical functions independent on the rel- 
ative velocity w, that corresponds to the vanishing coefficient Q = 0. This 
model describes a two-fluid system where the components are in full thermal 
isolation. Of course, each component can be a strongly-interacting matter 
rather than ideal gas. 

Indeed, a two-fluid model with the heat exchange between the components 
pertains better to superfluid systems because the "cold" and "warm" are 
not independent, and interference between them is evident. On the other 
hand, a two- or multi-component nuclear matter can be also described as 
an ideal two-component fluid if no macroscopic relative flow is expected. 
It can be a proton-neutron or quark-gluon system in equilibrium where only 
infinitesimal relative motion between the components is admitted in the form 
of spin-isospin sound. 

The the essence of two-fluid model is expressed in the decomposition of 
statistical sum (I55|) - (l57j) that results in additive formulas for the pressure 
(I59p and energy density (1601) . The statistical sum of the first component is 
calculated by standard formula fl5T|) . while the statistical sum of the second 
component is determined by formula fl52l) in the reference frame co-moving 
the first component. As a result, the relative velocity between the compo- 
nents w appears in the thermodynamical functions. 

When the pressure and energy density of interacting system is given f[6"2"]) - 
fl6"5|) . the two-fluid decomposition fl5T?l)- fl6ll is completed immediately accord- 
ing to formulas and (|6T))-(|68|). 

Both effects of coupling between fluid constituents and interaction be- 
tween particles and taken into account. The interaction between particles 
is responsible for non-ideality of the equation of state. It yields the poten- 
tial energy in addition to the kinetic energy of non-interacting system. In 
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the frames of mean-field approach it sets the effective mass M* and effective 
chemical potential /i* that appear in formulas (j5"2"l) - (155)1 . However, the hy- 
drodynamic non-ideality of this fluid system is not directly associated with 
the interaction between particles. 

Without loss of generality, the two-fluid decomposition fl59|) - fl60|) and fl65|) - 
f lB"8"j) can be also applied to superfluid systems, where the potential energy of 
interaction includes additional anomalous fields responsible for the conden- 
sation pressure and energy of the superfluid ground state [16]. For weakly 
interacting medium, like superfluid helium, the second (or "warm") compo- 
nent is often considered as an ideal gas of thermal excitations. However, such 
concept cannot be applied to superfluid nuclear matter because the contribu- 
tion of meson fields is not small, and exact two-fluid decomposition ( 16"5]) - (I6"81 
should be developed. It is the subject for further research. 

The authors are grateful to Erwin Schmidt for discussions. 
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